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ABSTRACT 
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Flutter  Analysis  of  a Cascade  of  Staggered  Blades  in 

Subsonic  Flow.  (December  1977) 

Louis  Kronenberger,  Jr.,  B.S.,  Texas  A&M  University 

Co-Chairmen  of  Advisory  Committee:  Dr.  Balusu  M.  Rao 

Dr.  Leland  A.  Carlson 


The  purpose  of  this  report  is  to  utilize  a numerical  lifting 
surface  theory  developed  by  Rao  and  Jones  to  predict  the  unsteady  air- 
loads on  an  infinite  cascade  of  staggered  blades  in  subsonic  compress- 
ible flow.  An  investigation  is  conducted  to  determine  the  effect  on  the 
unsteady  airloads  when  parameters  such  as  reduced  frequency,  interblade 
stagger  angle,  interblade  spacing,  and  interblade  phase  lag  are  varied 
over  a specific  range  of  values.  Once  the  unsteady  airloads  have  been 
determined,  they  are  used  to  perform  an  aeroelastic  analysis  of  the 
staggered  cascade  for  a single  degree  of  freedom  in  torsion  and  a two 
degree  of  freedom  system  in  bending  and  torsion. 

Results  of  the  single  degree  of  freedom  analysis  yields  flutter 
boundaries.-Iljese  are  compared  to  results  obtained  by  Whitehead  who 
used  a different  technique  for  calculating  the  unsteady  airloads  on  a 
finite  cascade. -^A  new  general  flutter  program  is  developed  for  the  two 
degree  of  freedom  system.^lh^airloads  are  used  as  forcing  functions  in 
the  resulting  two  Lagrangean  equ^ions  of  motion  representing  the  bend- 
ing and  torsional  degrees  of  freedomAfhe  iterative  procedure  of  the 


flutter  program  yields  the  flutter  frequency  and  speed  of  the  cascade 
reference  airfoil  as  a function  of  the  cascade  parameters. 
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NOMENCLATURE 

= speed  of  sound,  ft/sec. 
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= influence  coefficient 
= half-width  of  airfoil  segment 
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roots  of  the  quadratic  equation  for  the  real  part  of  +he 
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exponential  series  term 

3Sq/3X 

time,  sec. 

kinetic  energy  per  unit  span  of  airfoil 

3ct)/3x  = perturbation  velocity  component  along  x axis 

local  velocity,  ft/sec.;  blade  potential  energy 

3<}i/3z  = perturbation  velocity  component  along  z axis 
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= percent  blade  chord;  blade  axis  position 
= 

= cascade  interblade  stagger  angle,  deg. 

= <S/2v 

= (lo  + e)  = 
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= cascade  stagger  angle  as  defined  by  Whitehead  ; 
displacement  for  flapped  blade  in  transformation  of  axis 

= air  mass  density;  slug/ft^. 

= cascade  interblade  phase  lag 

= velocity  potential,  ft^/sec.;  phase  lag  as  defined  by 
Whi tehead 

= velocity  potential  at  a point  p 


= transformed  velocity  potential 
= pi/U  = reduced  frequency 


= critical  frequency 
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- (K^/M)^^  = blade  bending  natural  frequency 
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NOMENCLATURE  (CONT'D.) 

0)^  = (K^/I)^  = blade  torsional  natural  frequency 

V = gradient  operator  in  Cartesian  coordinates 

Superscripts 

j 

! ' = indicates  non-dimensional  ization  with  respect  to  the 
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. = d(  )/dt  = first  derivative  with  respect  to  time 

= d2(  )/dt2  = second  derivative  with  respect  to  time 

Subscripts 
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te  = blade  trailing  edge 


1 


INTRODUCTION 

Many  researchers  have  devoted  considerable  effort  towards  predicting 
the  flow  characteristics  through  the  multiple  stages  of  blades  that  exist 
in  axial  flow  compressors.  The  research  efforts  have  been  concentrated  on 
obtaining  the  aerodynamic  loads  utilized  in  the  design  of  more  efficient 
blades.  In  addition  to  the  need  for  information  about  blade  aerodynamic 
loading,  flutter  characteristics  also  need  to  be  determined. 

Assuming  two-dimensionality,  the  flow  through  a single  row  of  blades 
is  mathematically  equivalent  to  the  flow  through  a staggered  cascade  of 
infinitely  many  airfoils.  A staggered  cascade  is  shown  in  Fig.  1,  where 
X is  interblade  stagger  angle,  a is  blade  semi-chord  length,  U is  free- 

stream  velocity,  and  s is  interblade  spacing. 

1 2 

Early  researchers  like  Whitehead  , Kemp  and  Sears  , and  Schorr  and 

3 

Reddy  , assumed  incompressible  flow  and  developed  aerodynamic  theories 
and  computational  schemes  for  predicting  the  unsteady  airloads  on  the 

blades  for  oscillatory  freestream  flow  and/or  oscillating  blades. 

2 

An  example  is  the  work  of  Kemp  and  Sears  who  studied  the  problem 
of  the  unsteady  lift  generated  on  a reference  airfoil  of  a cascade.  In 
their  approach,  they  used  an  oscillatory  freestream  flow.  Their  study 
considered  the  steady  interaction  between  blades  but  neglected  the  un- 
steady interaction  and  therefore,  the  effects  of  cascade  spacing.  Their 
approach  was  to  express  the  unsteady  lift  as  a function  o-''  the  design 
parameters,  such  as  the  ratio  of  the  airfoil  chord  and  the  disturbance 
wavelength,  thus  enabling  a designer  to  optimize  the  performance  of  a 

The  format  of  this  thesis  follows  that  of  the  AIAA  Journal. 
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turbomachine  design  instead  of  analyzing  one  particular  blade  arrange- 
ment. 

3 

In  another  study,  Schorr  and  Reddy  treated  the  flow  through  a 
staggered  cascade  of  airfoils  in  which  the  effect  of  unsteady  upstream 
disturbances  were  included  as  in  the  case  of  unsteady  or  distorted  inlet 
flow  conditions  in  an  axial  flow  compressor.  Their  problem  was  also 
formulated  under  the  assumption  of  incompressible  potential  flow  and 
numerical  results  were  obtained  for  oscillatory  flow  using  an  approximate 
solution  developed  from  the  integral  equations  involved.  In  addition, 
their  solution  yielded  unsteady  lift  coefficients  for  the  airfoils  as  a 
function  of  the  frequency  of  the  oscillations  and  for  different  values 
of  stagger  angle  and  solidity  of  the  cascade. 

4 

In  an  independent  study,  Jones  and  Moore  studied  the  incompressible 
flow  about  a cascade  of  oscillating  airfoils  at  zero  mean  incidence.  In 
order  to  obtain  a solution,  they  utilized  a unique  numerical  lifting 
surface  technique  which  differs  from  all  other  methods  in  that  it  makes 
use  of  the  velocity  potential  instead  of  acceleration  potential  doublet 
distributions.  A major  advantage  of  this  method  is  that  it  leads  to  a 

5 

simpler  set  of  linear  equations.  Rao  and  Jones  later  applied  this  tech- 
nique to  the  oscillatory  flow  about  an  airfoil  of  a staggered  cascade. 
Airload  results  were  obtained  for  several  values  of  high  frequency, 
interblade  stagger  angle,  and  interblade  spacing  and  showed  excellent 

3 

agreement  with  the  results  of  Schorr  and  Reddy  . Results  were  also 
obtained  for  one  combination  of  interblade  stagger  angle  and  interblade 

spacing  at  several  interblade  phase  lag  angles  and  frequencies. 

3 6 

In  a technique  similar  to  that  used  by  Schorr  and  Reddy  , Fleeter 


used  oscillating  inflow  and  considered  the  effects  of  compressibility  on 
both  the  fluctuating  lift  and  fluctuating  moment  coefficients  for  cas- 
caded airfoils  having  an  upstream  non-uniformity.  A solution  was  obtained 
for  the  time-dependent,  two-dimensional,  partial  differential  equations 
which  describe  the  perturbation  velocity  potential  through  an  application 
of  Fourier  transform  theory.  The  resulting  integral  solution  equation 
was  evaluated  numerically  by  a matrix  inversion  technique.  The  fluct- 
uating lift  and  moment  coefficient  variations  were  computed  and  repre- 
sented as  a function  of  Mach  number,  cascade  solidity,  cascade  stagger 
angle,  interblade  phase  lag,  and  reduced  frequency. 

7 

Jones  and  Moore  extended  the  velocity  potential  formulation  to 
oscillating  two-dimensional  airfoils  in  compressible  flow.  In  their 
numerical  method  they  replaced  the  series  of  Hankel  functions  by  a rapid- 
ly convergent  exponential  series.  They  studied  the  effects  of  varying 
airfoil  spacing,  frequency,  Mach  number,  and  phase  difference  between 
adjacent  blades.  Variations  in  the  aerodynamic  damping  can  become  zero 
but  never  negative  at  certain  discrete  frequencies.  This  is  a desirable 
characteristic  with  respect  to  flutter  due  to  bending.  The  results  also 
indicated  that  the  pitching  moment  aerodynamic  damping  relative  to  the 
blade  quarter-chord  axis,  while  also  being  zero  at  the  critical  fre- 
quencies, could  be  negative  at  the  higher  Mach  numbers  over  a wide  range 
of  frequencies  of  interest  in  flutter  analysis.  This  is  an  undesirable 

characteristic  from  the  standpoint  that  it  increases  the  area  of  insta- 

8 

bility  for  torsional  flutter.  In  a recent  paper,  Rao  and  Jones  utilized 
the  theory  developed  in  Ref.  7 to  determine  the  airload  and  moment 
coefficients  on  a typical  airfoil  of  a staggered  cascade  of  airfoils  in 


subsonic  flow.  Circumferential  distortion  due  to  inflow  conditions  was 
expressed  as  an  interblade  phase  lag.  The  separate  conditions  for  oscil- 
latory inflow  and  oscillating  blades  were  considered.  Results  were 
obtained  for  several  values  of  frequency,  Mach  number,  interblade  spacing, 
interblade  stagger  angle,  and  interblade  phase  lag.  The  oscillatory  in- 
flow results  compared  well  with  those  of  Fleeter. 

9 

In  his  study  utilizing  compressible  flow,  Whitehead  presented 
calculations  for  the  torsional  flutter  of  a cascade  of  unstalled  blades 
at  zero  mean  deflection  and  subsonic  Mach  numbers.  Whitehead  found  that 
the  effect  of  increasing  Mach  number  was  favorable  and  tended  to  sup- 
press the  flutter  that  was  predicted  by  incompressible  theory. 

In  this  report,  a numerical  lifting  surface  theory  developed  by 

8 

Rao  and  Jones  is  used  to  predict  the  unsteady  airloads  on  an  infinite 
cascade  of  staggered  blades  in  subsonic  compressible  flow.  The  effect  of 
the  geometric  and  flow  parameters  on  the  airloads  is  investigated  by 
varying  them  over  a practical  range  of  values. 

Additionally,  an  investigation  is  conducted  for  a single  degree  of 
freedom  system  in  torsion.  The  effect  of  the  flow  and  geometric  para- 
meters is  evaluated  in  establishing  flutter  boundaries  and  these  results 

9 

are  compared  with  those  of  Whitehead  who  used  a completely  different 
computational  procedure  for  calculating  the  aerodynamic  derivatives. 

Finally,  a general  flutter  program  is  developed  for  a two  degree  of 
freedom  (bending- torsion) , staggered  cascade  in  subsonic  flow. The  geo- 
metric portions  of  the  equations  of  the  blade  motion  are  derived  using 

10 

Lagrange's  equation  of  motion  . The  unsteady  airloads  are  used  as 
forcing  functions  in  the  resulting  two  Lagrangean  equations  of  motion 


6 


representing  the  bending  and  torsional  degrees  of  freedom.  By  utilizing 
an  iterative  procedure  which  permits  frequency  variation,  the  flutter 
frequency  and  flutter  speed  of  the  reference  airfoil  are  obtained  as  a 
function  of  the  cascade  parameters. 
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AERODYNAMIC  THEORY 


General 


The  governing  equation  for  the  unsteady,  conipressihle,  two- 

dimensional  flow  of  an  iscntropic,  inviscid,  irrotational  fluid  is 
1 1 

given  in  terms  of  its  velocity  potential  by. 


\'U 


1 


= 0 


(1) 


where  t is  the  perturbation  velocity  and 

0 = {U+u)i  + wk  . (2) 

Freestream  velocity  is  given  as  U.  The  respective  perturbation  velocity 

compenents  along  the  x and  z axes  are,  li(=||-)  and 

Assuming  that  u and  w are  very  small  compared  to  U,  Fqs.  (1)  and 
(2)  are  combined  to  yield. 


ax-^  ^ 3z-^  - a2  a#  ^ T alcat  ’ 

where  M is  the  freestream  Mach  number  and  a is  the  speed  of  sound. 

7 

Jones  used  a non-dimensional  coordinate  transformation  such  that, 

X = ii.  Z = fi-  T = — 
t £ £. 

and 

♦(x.z.t)  = U»(X,Z)e'(^X-T) 

where  £ is  a reference  length  for  blade  semi -chord  and  & = (1-M^)'*. 
Also, 


(3) 


(4) 


(5) 


= ei 
u ’ 


c - 


, and  t = v(X ,Z)  . 


(6) 


When  Eq.  (3)  is  conbined  with  Eq.  (5),  it  reduces  to  a two-dimensional 
Helmholtz  equation  for  the  perturbation  velocity  in  the  transformed  co- 
ordinate system, 


V2(j.  + = 0 


(7) 


, Mo) 

where  = — 


For  a flow  problem,  the  boundary  conditions  are  usually  prescribed. 

For  an  isolated  airfoil  in  compressible,  unsteady  flow,  the  solution  of 

12 

Eq.  (7)  can  be  derived  by  the  application  of  Green's  Theorem  . A rela- 
tion for  the  velocity  potential  at  a given  point,  4.^,  is  given  in  terms 
of  velocity  potential  distribution  over  the  lifting  surface  and  its 
wake.  Treating  the  lifting  surface  as  a thin  airfoil,  the  discontinuity 
in  the  velocity  potential  between  the  upper  and  lower  surfaces  is  ex- 
pressed as  a doublet  distribution,  K(=  “^^pper  ~ ‘^lower^’  solu- 

tion for  an  isolated  airfoil  is  given  in  Ref.  7.  The  relation  between 
the  downwash  at  any  point  p on  a thin  reference  blade  on  a staggered 
cascade  in  subsonic  flow  and  the  modified  doublet  distribution  K(X),  is 
given  as. 


2t.W  = - K(X)^-^_  S (X  -X,Z  ,D,H,o)  dX  (8) 

p J ' ' ? o'  p ’ p’  ’ ’ ' 

-1 

where , 

CD 

H^^^|^r[(X^-X+mD)2  + (mH-Zp)2]’''=j  (9) 

m=-“ 

and  D ^ ^ ''s  the  interblade  phase  lag.  The  blades  of 

the  cascade  are  numbered  m,  with  m = 0 as  the  reference  airfoil.  Since 


9 


3^5  3^5 

and  in  the  limit  as  -+  0,  Eq.  (10)  may  be  rewritten  in  the  form. 


(10) 


2r,W(Xp)  = 


where , 


S,  = 


-1 


/3S 

«x|ax-  ’ « 


- + K?S 


dX 


(11) 


3S, 


1 3X 


(X  -X+mD)Hj^V 


0 , jLk\  pini(o+cD) 


(X^-X+mD)24m2H2] 


m=-“' 


[(X  -X+mD)2+m2H2j’= 


.(12) 


The  series  involving  Hankel  functions  (S^,S^)  in  Eqs.  (9)  and  (12)  have 

8 

very  poor  convergence  characteristics  . Therefore,  these  are  replaced 
by  an  exponential  series  as  shown  in  Refs.  13  and  14.  However,  it  is 
important  to  understand  that  this  transformation  from  Hankel  function 
series  to  exponential  series  is  valid  only  for  an  infinite  cascade. 
Hence,  it  cannot  be  applied  to  a finite  cascade.  The  convergence  of  the 
exponential  series  is  so  good,  that  the  required  computational  time  is 
less  than  for  a cascade  when  compared  to  a two-dimensional  isolated 
airfoil  in  subsonic  flow  where  it  is  required  to  use  Hankel  functions. 
The  transformed  relations  as  given  in  Ref.  7 are. 


A -2^a(m)|X  -Xl/S 

s -_ire p 

[(6-m)2-u2j2 


(13) 


ni=-o 


and , 


^1  = ^ 


-2Tia(m)|X  -X|/S 

a(m)e 


m=-° 


[(6-m)2-u2j’^ 


for  X > or  < X. 
P 


(14) 


where , 


10 


I 

i 


I 


a(m)  = [(6-m)2-u2]*^  ^ ± i(6-m)  ^ , for  Xp  > or  < X,  (15) 

and, 

6 = 2^  , S = {D2+h2)^"  . (16) 

The  exponential  form  of  the  series  and  not  only  converge 

rapidly,  but  also  provides  directly  the  values  for  critical  frequencies 
as  shown  in  Ref.  7.  The  Eqs.  (13)  and  (14)  will  diverge  whenever  one  of 
the  denominators  vanishes  and  no  solution  of  Eq.  (11)  would  be  possible. 
The  critical  values  of  the  parameter  u for  which  the  analysis  fails  are 
given  by, 

p=5,  1±6,  2±6,  etc. 

This  phenomenon  corresponds  to  a resonance  condition  at  critical  fre- 
quencies which  constitute  an  infinite  set  of  values  of  a parameter 
depending  on  flow  and  configuration  characteristics  at  which  the  aero- 
dynamic function  becomes  infinite  everywhere.  Resonance  conditions,  as 
shown  in  App.  A,  are  functions  of  Mach  number,  frequency,  interblade 
spacing,  phase  lag,  stagger  angle,  and  acoustic  velocity.  They  repre- 
sent the  condition  at  which  self-induced  aerodynamic  forces  are  zero 
and  the  blades  act  effectively  as  if  they  were  in  a vacuum. 

Boundary  Conditions 

The  downwash  w(=w'e^'“'^)  can  be  expressed  in  terms  of  il>(X,Z)  by 
using  Eq.  (5), 

w = If  = If  . (17) 


Downwash  can  now  be  non-dimensional i zed  to  give. 


11 


' -dl  ~ bu 


Rao  and  Jones  studied  the  effects  of  unsteady  airloads  on  a cas- 
cade of  staggered  blades  in  subsonic  flow  for  both  oscillatory  inflow 
and  oscillating  blades.  For  the  oscillatory  inflow  condition,  Rao  and 
Jones  assumed  a sinusoidal  gust  loading  such  that  the  downwash  boundary 
conditions  on  the  airfoil  were. 


w.  = u e ^ 


where  i = 1,2,....N.  For  the  oscillating  blades  with  flapping  and 
pitching  motions  about  the  mid-chord  position  as  shown  in  Fig.  2,  the 
downwash  boundary  condition  is  defined  as. 


I p I i - 

w^  = U i(i)Z  + { 1 +i  )ci 


where  z and  a are  the  amplitudes  in  flapping  and  pitching  respec- 


tively. Since  periodic  motions  were  assumed,  z and  a are  defined  as. 


ipt  _ iwT 
z = 2e'^  = ze 


’ ipt  ' iwT 
a=ae  =ae 


(21a) 


(21b) 


The  K^'s  are  normally  complex  and  depend  on  z and  a for  any 
particular  values  of  u and  M.  A typical  will  have  the  form, 

I I 

K=az+ba  (: 

n n n ' 

where  a^  and  are  complex  quantities  and  depend  on  frequency,  Mach 

number,  interblade  spacing,  interblade  stagger  angle,  and  interblade 
phase  lag. 
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Nuiiiericdl  Procedure 


The  numerical  lifting  surface  method  as  developed  in  Ref.  7 is 
utilized.  The  wake  boundary  condition  is  given  by. 


K(X)  = K.  e""" 
te 


(23) 


where  is  the  doublet  strength  at  the  trailing  edge  (X=l).  Wiien  fqs. 


(11)  and  (23)  are  combined,  they  can  be  expressed  in  the  form, 

^ 1 


2ttW 


(24) 


where , 


-1 


It 


ivSot  - v^(1-M2)p 


P = 


^g-2Tia(m)(l-X^)/S 


m=-“i 


(f-M)2-p2T' 


(25) 

(26) 


|2TTa(m)+i  vS 

and  a(m)  = j^(6-m)^-p^J  ^ - i(6-m)  ^ . The  symbols  and  S.|^  represent 

and  S.J  respectively  at  the  trailing  edge. 

The  airfoil  is  divided  into  N equal  strips  and  K is  assumed  to  be 


, a constant  over  the  nth  strip,  where  n = 1,2 N.  Let  2B  be  the 


width  of  each  strip  and  X^  denote  the  center  of  the  nth  strip.  In  Eg. 


(24),  Xp  and  are  replaced  by  X.  and  respectively,  where  i refers 


to  the  ith  strip.  Hence,  Eg.  (24)  will  then  be  given  approximately  by, 
N 


27iW.  =)  K Ts/X.-X  -B)  - S.(X.-X  +B)  + 2*,=’BS  (X.-X  )] 
1 / n L 1 1 n ' 1 ' 1 n ' o ' i n ' J 


n=l 
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+ K^^I(X^)  , where  i = 1,2 N.  (27) 

When  = X^,  is  evaluated  by  integrating  the  exponential  form  of 
over  the  interval  X^-B  < X < X^+B.  For  otiier  intervals,  X^-  f X^  and 
can  either  be  found  by  integration  or  be  taken  as  the  mean  value  of 
over  the  interval.  K,.  can  be  expressed  as  a function  of  K by  usina  the 
wake  condition.  This  relation  is  given  by. 


K 


te 


K 

n 


(2ivB+e 


(28) 


For  a given  geometry  and  flow  conditions,  the  numerical  terms  in 
Eq.  (27)  can  be  evaluated  by  using  Eqs.  (13),  (1^),  (25),  and  (26).  For 
a known  set  of  WVs,  Eq.  (27)  reduces  to  a system  of  N equations  with  N 

unknowns,  K.| , K^, K^,  when  it  is  combined  with  Fq.  (23).  It  is 

assumed  that  is  known  over  the  airfoil.  In  Eq.  (27),  this  (^presents 
a set  of  linear  algebraic  equations  given  by, 

27i{W.}  = fA]{K}  . (29) 


Therefore,  knowing  the  values  of  {W.. } and  [A]  allows  for  the  determina- 


tion of  the  doublet  distribution,  {K}  . 

Aerodynamic  Derivatives 
1 s 

Euler's  equation  of  motion  is  given  as, 

3U  ^ 1^90  _ J_  9_P 
'9t  ^9X  ' p '9X  ■ 


(30) 


In  terms  of  the  upper  and  lower  velocity  potential  and  pressure  on  the 
airfoil,  Eq.  (30)  can  be  shown  to  be, 
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9k  , I (3  k 1 / r.  rN  \ 

where  k(x)  = • Since  lift  per  unit  chord  per  unit  span  is 

i(x)  = (PJ,-P^J).  Eq.  (31)  becomes. 


Now  from  Eq . ( 5 ) , 


k(x)  = 


k(x)  = , 


(33a) 


(33b) 


where  K(X)  = (‘!’^j"‘J’j^)  • Equations  (32)  and  (33)  are  combined  to  yield, 

t(X)  = pU2[ivK(X)  + i^Xjjgi(cX+aiT)  (3^) 

where  v = w/6^.  Equation  (34)  is  valid  on  the  airfoil  surface, 

-1  i X < 1.  In  the  wake  region,  no  pressure  discontinuities  are  allowed 


to  exist  and  hence. 


vK(X)  + = 0 


must  be  satisfied  when  X ;;  1 . 

When  values  of  have  been  obtained,  the  local  lift,  i,  at  a 
point  X is  given  by  Eq.  (34).  If  K is  substituted  for  Ke  then 


^ = (ia,l<  + |^)e’“'^ 


(36a) 


~^=  (iu,K  + |^)e^P^ 


(36b) 


The  lift  L(=L  e^P^)  and  the  nose-up  pitching  moment  about  the  mid- 
chord axis  M(=y  e^P^)  are  given  by. 


\ 

L 
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r 


and , 


- K.  + i u 

pU2£  ^ 


K dX  = C,  z + C a 

iiZ  9.a 


(37) 


■1 


pU^i 


= K.  + 
2 te 


r 


K dX  - ii 


1 

XK  dX 


- C z + C Cl 

mz  iTici 


(38) 


where  C.  , C„  , C , and  C are  the  aerodynamic  derivatives.  These  are 
£.2  ta  inz  nia  ■' 

usually  complex  and  depend  on  the  geometry  and  flow  characteristics. 


Tquations  of  Motion 


Flutter  can  be  defined  as  the  dynamic  instability  of  an  elastic 
body  in  an  airstream.  A dynamic  system  is  assumed  to  flutter  at  a speed 
where  the  net  damping  forces  are  equal  to  zero.  There  exists  a frequency 
at  which  flutter  occurs.  This  is  termed  the  flutter  frequency.  Consider 
the  spring  supported  airfoil  given  in  Fig.  3 where  and  are  stiff- 
ness in  bending  and  torsion,  respectively.  The  blade  segment  is  permit- 
ted freedom  to  execute  small  periodic  vertical  displacements, 

iz{=lz  e^'*^^),  and  angular  displacements,  a{=a  e^^^),  where  U is  free- 
stream  velocity  and  an  oscillation  frequency  of  p rad/sec.  The  cascade 
parameters  associated  with  a staggered  row  of  thin  oscillating  blades 
are  presented  in  Fig.  1. 


For  the  reference  airfoil,  the  expression  for  the  kinetic  energy 
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where  z and  a are  nondimensional  coordinates.  By  integrating  Eq.  (39), 
the  kinetic  energy  per  unit  span  of  the  airfoil  is  given  by, 


where , 


T = ^ (tz  + ra)^dm 

J 

1 1 
T = 2 Mt^Z^  + SJZa  + 2 la^ 


M = dm  = blade  mass  per  unit  span. 

S =J^r  dm  = blade  static  moment  about  elastic  axis. 

I = ^ r^  dm  = blade  mass  moment  of  inertia  about  elastic  axis. 


The  potential  energy  expression  for  the  blade  is  given  by, 

U = t K £,2z2  + } K . ( 

2 z Z Cl 

The  bending  and  torsional  stiffness  coefficients  are  related  to  the 
respective  natural  frequencies,  and  by  the  relations. 


(43a) 


K = . 

a a 


(43b) 


Flutter  is  a self-exciting  aeroelastic  phenomenon  involving  the  inter- 
action between  the  inertial,  elastic,  and  aerodynamic  forces.  The 
structural  damping  of  the  blade  can  be  neglected  compared  to  the  aero- 
dynamic damping  terms  associated  with  the  forcing  function.  Neglecting 
the  structural  damping,  Lagrange's  equations  of  inotion  can  be  expressed 


JT 

dt 


. = F. 


s 
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where  F|^  and  q|^  are  generalized  forces  and  coordinates.  The  aerodynamic 

forcing  functions  are  obtained  from  the  two-dimensional  unsteady  air- 
loads prediction  method  and  are  given  in  the  form, 

(45a) 


(45b) 


where  , and  are  the  complex  aerodynamic  derivatives 

referred  to  the  elastic  axis  and  are  given  in  real  and  imaginary  terms 
as, 


C ~C  .,n  + iC  , 

zz  t’R  tzi 


C ~Cn+iC„T 
Za  tfiR  £al 


^mz  " ^mzR  ^^mzl 


C “ C n i c T 
HM  ITKtK  mxi 


(46a) 

(46b) 

(46c) 

(46d) 


Now  Fqs.  (41),  (42),  and  (45)  are  combined  with  Eq . (44)  to  obtain, 


mz 


+ C„  a) 

ZZ  Za 

(47a) 

Z t C a)  . 
nu 

(47b) 

The  solutions  of  Eqs.  (47)  are  assumed  to  have  periodic  motion  given  by, 

(48a) 


' ipt 
z = z e 


' ipt 
a = a e 


(48b) 


Substitution  of  Eqs.  (48)  into  Eqs.  (47)  and  combining  terms, 


(tK^  - Mtp2  + pU2jiC^^)z  + - Sp2)a  = 0 


+ S£.p^)z  + (K^  - _ ip2)„  = 0 


(49a) 

(49b) 


A 
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Equations  (49)  are  a set  of  two  simultaneous,  linear,  homogeneous, 

I I 

algebraic  equations  in  2 and  u . For  this  system  to  have  a non-trivial 

I I 

solution,  the  coefficient  determinant  of  z and  a must  be  equal  to 


zero. 


(tK^  - Mep2  + pU2£C^^) 


(pU2cCj,^  - Sp  ) 


-(C  pU2j,2  + S)ip2) 
' mz 


- Ip2) 
ma  ^ 


= 0 


(50) 


- »• 

This  determinant  is  referred  to  as  the  flutter  determinant  and  the 
solution  gives  the  flutter  frequency.  However,  a direct  solution  cannot 
be  obtained  since  the  aerodynamic  derivatives  are  functions  of  w(  = ps./U). 
The  aerodynamic  derivatives  are  complex  quantities  and  thus  the  flutter 
determinant  can  be  expressed  in  two  parts,  real  and  imaginary.  After 
expanding  the  flutter  determinant  and  substituting  appropriate  values  of 
Eqs.  (46),  the  real  part  becomes, 

[(HU  - iH)  - 4 


p2£7 


If-(C„  nC  D - c,  ,C  - + C„  .C  , - C„  dC  d)]  = 0 
' 07f?  rrVfR  CtT  mr»  T jloi  I IHZ  I OrvW  rn7W'J 


'nzR  rmR 

and  the  imaginary  part  becomes. 


'JiaR''mzR' 


(51) 


^ ^^^taR^mzI  ^JlciI^mzR  ‘ ^tzR^nviI  ‘ ^CzI^maR^l  ""  ° 

where  P = — .If  the  real  and  imaginary  parts  are  set  equal  to  zero, 
p2 
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li 

ji 

I' 

I, 

! 

I 
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solutions  PI  and  P2  can  be  obtained  for  the  real  part  and  P3  for  the 
imaginary  part.  For  an  assumed  value  of  w,  if  one  of  tlie  solutions,  PI 
or  P2,  of  the  real  part  is  equal  to  the  solution  P3  of  the  imaginary 
part,  then  this  w corresponds  to  the  flutter  frequency.  The  main  objec- 
tive is  to  find  the  flutter  speed.  However,  the  aerodynamic  derivatives 
are  functions  of  the  Mach  number  and  the  reduced  frequency.  The  flutter 
problem  can  only  be  solved  after  the  aerodynamic  derivatives  are  eval- 
uated. It  is  necessary  to  assume  a Mach  number  and  a reduced  frequency 
and  test  whether  flutter  occurs  at  these  values.  If  the  test  results  are 
negative,  then  it  is  necessary  to  iterate  on  reduced  frequency  until  a 
flutter  speed  is  obtained  for  the  assumed  Mach  number.  If  the  two  speeds 
are  not  equal,  then  it  is  necessary  to  iterate  on  Mach  number  until  the 
flutter  Mach  number  is  equal  to  the  assumed  Mach  number.  A computational 
example  is  illustrated  in  App.  B. 

The  reference  airfoil  as  illustrated  in  Fig.  3,  with  its  center  at 
the  origin  of  coordinates,  is  assumed  to  be  describing  plunging  and 
pitching  oscillations  of  frequency  p rad/sec.,  defined  by  z and  a.  For  a 
single  degree  of  freedom  of  flapping  reference  blade,  the  equation  of 
motion  is  given  by. 


Mz  + K z = L z 
z z 


where  z is  defined  for  periodic  motion  as. 


' ipt 
z = z e 


3Z 


ipt 

^ = ipz  . 


is  a complex  quantity  with  real  and  imaginary  components. 


(53) 


(54) 

(55) 
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L = L „ + iL  , 
2 zR  zl 


(56) 


Now  Eqs.  (55)  and  (56)  will  give. 


I z = zL  r,  + izL  , = zL  „ + - L , 
z zR  zl  zR  p zl 

and  Eq.  (53)  now  becomes. 


(57) 


* "'z  ■ ">  • '581 

when  like  terms  of  Eq.  (58)  are  compared  to, 

Mz  + C^z  + Kz  = 0 . (59) 

The  damping  constant,  becomes  a function  of  the  imaginary  part  of 
the  lift  coefficient  in  ^lapping. 

L , 

C,  = - --  (60) 

L p 

In  a similar  manner,  if  a single  degree  of  freedom  in  pitching  is 
considered,  the  equation  of  motion  is  given  by. 


!«  + K a = M a 
a Cl 


(61) 


where , 


' i pt 
a = a e 


' . i pt 

a - a 1 pe  - 1 pu 


(62a) 

(62b) 


and  considering  the  fact  that  is  also  a complex  quantity  given  as. 


M = M D + iM  , . 
Cl  aK  al 


(63) 


It  can  be  shown  that, 


'8  - f * “(K,  ■ n,R>  = » 


(64) 
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RESULTS  AND  DISCUSSION 

Convergence  Studies  and  Variation  of  Parameters 

A convergence  study  was  conducted  utilizing  a program  developed  by 
16 

B.M.  Rao  . This  program  used  the  numerical  lifting  surface  technique 
to  derive  the  aerodynamic  derivatives  for  a cascade  of  staggered  blades. 
In  the  analysis,  the  downwash  induced  at  the  reference  airfoil  by  the 
doublet  distribution,  K(X),  over  each  oscillating  blade  of  the  cascade 
is  represented  by  a series  of  Hankel  functions  denoted  by  S^.  This 
series,  in  the  form  presented,  is  poorly  convergent.  However,  it  has 
been  shown  in  Refs.  13  and  14  that  it  can  be  replaced  by  a rapidly  con- 
vergent series  of  exponential  terms.  The  problem  was  evaluated  using 
the  alternative  exponential  form  of  and  its  derivative  form  of 
S^{=3S^/8X). 

An  influencing  parameter  in  the  convergence  of  the  exponential 
function  is  the  number  of  terms  that  are  included  in  the  series  before 
it  converges  to  an  acceptable  value.  While  the  value  of  the  exponential 
series  is  dependent  on  the  sum  of  the  number  of  terms  considered,  a 

_5 

change  of  1x10  was  established  as  a constraint  criterion  between 
two  consecutive  terms  for  convergence  or  until  a prescribed  limit  on 
the  number  of  terms  included  in  the  series  was  reached.  In  the  program, 
NNN(=2L+1,  where  L is  any  integer  from  1 to  “)  was  the  control  para- 
meter that  was  varied  in  order  to  determine  the  maximum  limit  of  terms 
that  could  be  summed  in  the  series.  Naturally,  the  greater  number  of 
terms  would  tend  to  converge  on  a more  accurate  value  for  the  series. 


r 


25 


However,  in  the  interest  of  program  efficiency,  a coii.promise  value  for 
NNN  had  to  be  determined  and  accepted.  Table  1 shows  typical  results 
obtained  for  various  values  of  NNN  while  varying  w from  0.10  to  0.40, 
dividing  the  reference  blade  into  16  segments,  with  an  interblade  stag- 
ger angle  of  45  degrees,  Mach  number  of  0.60,  interblade  phase  lag  of 
90  degrees,  and  a blade  to  semi-chord  space  ratio  of  2.0.  As  is  shown  in 
the  table,  the  influence  of  NNN  on  the  aerodynamic  derivatives  when  it 
is  greater  than  100  becomes  insignificant.  Therefore,  a value  of  NNN 
equal  to  100  was  chosen  for  computational  purposes. 

An  additional  program  parameter  that  was  considered  in  the  con- 
vergence study  was  that  of  the  number  of  segments  on  a typical  airfoil. 
Consider  the  two-dimensional  airfoil  as  shown  in  Fig.  4.  The  blade  is 
divided  into  N segments  with  the  coordinate  axis  located  at  the  mid- 
chord of  the  blade  and  a collocation  point,  X^.  , located  at  the  center  of 

each  segment  ( where  i = 1,2, ,N).  The  doublet  distribution,  K(X), 

over  each  segment  is  assumed  to  be  constant. 


X^  = Collocation  points 


Fig.  4 Segmented  Blade  for  N = 4 . 
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Table  1 tffect  of  Parameter  Variation  on 
i Aerodynamic  Derivatives 


50 

JINN_ 

100 

150 

0)  = 0.10 

C 

niz 

C 

mu 

0.2895  + 0.6259i 

6.4762  - 2.6460i 

0.2356  + 0. 28111 

2.9262  - 2.4431i 

0.2895  f 0. 62601 

6.4766  - 2. 64671 

0.2357  + 0.28111 

2.9263  - 2.44391 

0.2895  + 0.62601 

6.4766  - 2.64671 

0.2357  + 0.28111 

2.9263  - 2.44391 

W = 0.20 

0.5423  + 1.0312i 

0.5425  + 1 .03121 

0.5425  + 1.03121 

^J!.a 

5.5959  - 2.3289i 

5.5963  - 2. 32981 

5.5964  - 2.32991 

^mz 

0.4771  + 0.4295i 

0.4773  + 0.42951 

0.4773  + 0.42951 

2.3739  - 2.6314i 

2.3739  - 2.63221 

2.3739  - 2.63231 

a.  = 0.30 

0.6376  + 1.3361i 

0.6378  + 1 .33621 

0.6379  + 1 .33621 

^Za 

5.0769  - 1.6554i 

5.0774  - 1 .65531 

5.0774  - 1.65651 

Saz 

0.6689  + 0.5233i 

0.6691  + 0.52331 

0.6691  + 0.52331 

^ma 

2.0558  - 2.6647i 

2.0558  - 2.66561 

2.0558  - 2.66571 

(0  = 0.40 

^tz 

0.5794  + 1.6506i 

0.5798  + 1.65081 

0.5799  + 1.65081 

C 

4.9109  - 0.8886i 

4.9117  - 0.88961 

4.9118  - 0.89001 

c 

ITIZ 

0.8300  + 0.61971 

0.8303  + 0.61971 

0.8304  + 0.61971 

c 

1.9325  - 2. 70851 

1 .9325  - 2.70961 

1.9325  - 2.70981 

nva 
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The  nuinber  of  segments  chosen  directly  influences  the  value  of  the 
integral  evaluated  to  determine  the  downwash  at  each  collocation  point. 
Since  a larger  numt)er  of  segments  will  converge  to  an  increasingly  more 
accurate  value,  a study  was  conducted  to  determine  the  most  significant 
number  of  segments  along  the  blade  to  be  used.  While  a minimum  of  four 
segments  could  be  used,  the  degree  of  confidence  in  the  downwash  values 
increases  for  segments  greater  than  N = 10.  The  program  evaluates  the 
downwash  integral  using  a Simpson's  three-eighths  rule  scheme.  Addition- 
ally, the  first  and  last  half-segment  are  evaluated  by  using  the  trape- 
zoidal rule.  This  requires  multiples  of  N(  = 3L+1  , wiiere  1.  is  an  integer) 
segments  to  be  used.  Table  2 shows  the  typical  influence  that  varying 
N has  on  the  aerodynamic  derivatives.  While  f = 16  provides  reasonable 
accuracy,  N = 22  was  chosen  to  be  used  throughout  the  program  compu- 
tations. Theoretically,  the  greater  the  numl)Gr  of  segments  chosen,  the 
greater  the  accuracy  of  the  aerodynamic  derivatives.  However,  because 
of  the  way  the  solution  technique  is  employed,  it  is  possible  that 
numerical  instabilities  could  be  present  for  values  when  N becomes 
much  greater  than  twenty-five. 

For  flutter  studies,  the  aerodvnamic  derivatives  C and  C 

a Z ITla 

provide  meaningful  information  and  are  shown  for  several  sets  of  the 

geometric  and  flow  parameters.  Figures  5 thru  12  show  the  out-of-phase 

components,  (C.  , and  C of  C,  and  C which  represents  the  damp- 
szl  mal  ez  mu  f 

ing  terms  in  pure  translational  and  pitching  motions,  respectively. 

The  sign  convention  used  is  such  that  the  negative  values  of  C 

m:i 

represent  the  positive  damping  constant  and  the  positive  values  of 


Table  2 Effect  of  Blade  Segment  Variation  on 
Aerodynamic  Derivatives 


0.4527  + 0.74911 
7.8324  - 4.2956i 
0.3850  + 0.2S38i 
3.0517  - 4.03511 


(V  = 0.20  I 

Sz 
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mz 
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0.8705  + 1.15711 

6.4617  - 4.09921 

0.7579  4 0.32551 

1 .8926  - 4.30111 

111  = 0 . 30 

' S,z 

'mz 

C 

m.Ti 

1 .1271  4 1 .47661 

5.8947  - 3.56701 

1 .0689  4 0.26471 

1.1571  - 4.46911 

uj  = 0.40 

C.Z 

Snz 

C 

nia 

1.3239  4 1.86431 

5.9616  - 3.26261 

1 .3968  4 0.1  31  31 

0.5074  - 4.86491 

uj  = 0 - 50 

Cmz 

1.6653  4 2.34491 

6.3849  - 3.65121 

1.7879  - 0.21371 
Lo.5799  - 5.51621 

0.4520  -t  0. 74561  0.4518  + 0.74361 
7.8069  - 4.27661  7.7898  - 4. 26941 
0.3796  + 0.27351  0.3764  + 0.27361 
2.9508  - 3.98351  2.9015  - 3.95331 


0.8637  + 1.15041  0.8604  + 1.14601 

6.4427  - 4.04651  6.4265  - 4.02201 

0.7453  + 0.31151  0.7381  + 0.30531 

1.8251  - 4.24491  1 .7945  - 4.21091 

1 .1  1 33  + 1 .471  51  1.1062  + 1 .46721 

5.8968  - 3.49601  5.8398  - 3.46131 

1.0540  + 0.25051  1.0454  + 0.24511 

1.1128  - 4.43091  1 .0961  - 4 .40501 

1 ..3050  + 1 .86381  1.2948  + 1 .86131 

5.9857  - 3.18351  5.9889  - 3.14321 

1 .3843  + 0.1  1451  1.3768  + 0.10941 

0.4636  - 4.85201  0.4579  - 4.83831 

1.6441  -i  2.35101  1.6321  + 2.35221 

6.4326  - 3.57191  6.4488  - 3.52941 

1.7789  - 0.23931  1.7742  - 0.24631 

0.6328  - 5.52491  -0.6457  - 5.52441 
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Cj,^j  represent  positive  translational  stability. 

Figures  5 and  6 show  the  compressibility  effects  (M  = 0.3,  0.5, 

0.7,  and  0.8)  on  C and  C when  s = ?e,  A = 45^’’,  and  o = 180*^  as  co  is 
kZ 

varied  from  0.1  to  0.5.  The  absolute  values  of  are  shown  to  increase 

as  w increases  for  Mach  numbers  less  than  M =0.7  and  undergoes  rapid 

changes  for  M = 0.8.  For  the  parameters  considered  at  M = 0.8  in  Fig.  5, 

there  is  a rapid  change  for  values  between  w = 0.40  and  w = 0.51.  The 

lowest  critical  frequency  occurs  at  w = 0.51.  Figure  6 shows  a tendency 

toward  larger  changes  at  the  higher  M = 0.8  value.  As  the  Mach  number 

increases,  the  maximum  value  of  damping  also  increases.  This  indicates 

a stabilizing  influence  in  damping  as  Mach  number  increases. 

Figures  7 and  8 show  the  effect  of  interblade  spacing  (s  = 1.6£, 

2.0i,  and  2.4?,)  on  C„  , and  C ,,  respectively  when  M = 0.8,  A = 45°, 

tzi  niul  I j 

and  0 = 180°  as  u>  is  varied  from  0.1  to  0.5.  For  the  parameters  and  +he 
range  of  frequencies  considered,  the  critical  frequencies  are  for 
s = 1.6?,  = 0.64;  s = 2.0?,  = 0.51;  and  for  s = 2.4?,  = 0.42. 

Both  figures  indicate  that  as  the  interblade  spacing  increases,  the 
lowest  critical  frequency  decreases.  As  uj  increases,  the  damping 
decreases  rapidly  to  zero  at  u = and  then  increases  again  as  w is 

increased  without  becoming  negative.  Below  the  first  critical  frequency, 

the  damping  decreases  when  the  spacing  between  the  blades  is  increased. 

Figures  9 and  10  show  the  effect  of  blade  stagger  angle  (A  = 0°, 

45°,  and  60°)  on  C , and  C ,,  respectively  when  M = 0.8,  s = 2.0?, 

? z 1 mu  1 


and  0 = 180°,  as  w is  varied  from  0.1  to  0.5.  For  the  parameters  and 
the  range  of  frequency  values  considered,  the  critical  frequency  values 
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are  for  A = 0°,  lu^  = 1.18;  A = 45°,  = 0.51;  and  A = 60°,  = 0.44. 

As  tlie  stagger  angle  increases,  the  lowest  critical  frequency  decreases. 
All  the  coefficients  are  presented  with  respect  to  the  reference  axis  at 
mid-chord.  In  Ref.  7,  Jones  has  shown  that  it  is  possible  to  incite  pure 
pitching  oscillation  in  the  airfoils  in  cascade  when  A = 0°  and  0.73  < u> 

< 0.98  with  a reference  axis  forward  of  the  inid-ciiord.  However,  since 
the  values  given  in  Fig.  9 only  consider  the  practical  range  of  values 
for  0.1  < (D  < 0.5,  the  translational  damping  coefficient  remains 

positive  indicating  that  there  is  no  pure  translational  instability. 

Figures  11  and  12  show  the  effect  of  interblade  phase  lag  (o  = 0°, 
90°,  and  180  ) on  and  , respectively  when  M = 0.8,  s = 2.0)!., 

and  A = 45°  as  w is  varied  from  0.1  to  0.5.  For  the  parameters  and  the 
range  of  frequency  values  considered,  the  critical  frequency  for  o = 180° 
is  = 0.51.  The  translational  and  pitching  damping  constants  remain 

positive  throughout  the  range  of  frequency  values  considered.  However, 

as  a is  increased  with  increasing  the  damping  decreases  rapidly  for 

0 = 180°  until  it  equals  zero  at  w = u . 

c 

Single  Degree  of  Freedom 

A single  degree  of  freedom  analysis  for  torsion  was  performed  on 
a cascade  of  blades  using  the  lifting  surface  theory  developed  by  Jones 
and  Rao.  After  obtaining  the  results,  a comparison  was  made  to  the 

9 

results  obtained  by  Whitehead  who  employed  a technique  developed  by 

1 7 

Smith  to  obtain  the  aerodynamic  forces  and  mofnents.  Whitehead  dis- 
tinguished a region  of  flutter  that  he  termed  as  "sub-critical  flutter”. 
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(i  obtained  for  a space-to-chord  ratio  of  1.0,  a stagger  angle  of  45 

degrees,  a position  of  the  torsional  axis  at  58.8  percent  cnord,  and  an 
interblade  phase  angle  of  60  degrees.  The  imaginary  part  of  the  iriuiTient 
coefficient  is  plotted  against  the  reduced  frequency  parameter.  Also 

I shown  in  Fig.  13  is  a plot  of  the  imaginary  part  of  the  moment  coeffici- 

[ 

I ent  as  obtained  by  utilizing  the  Jones/Rao  technique.  As  can  be  seen, 

I when  given  the  same  conditions,  the  results  obtained  when  using  the 

^ Jones/Rao  lifting  surface  technique  agree  closely  with  the  results 

i obtained  by  Whitehead.  As  shown  previously,  the  imaginary  part  of  the 

moment  coefficient  is  capable  of  adding  energy  when  it  is  negative. 
Assuming  that  there  is  no  mechanical  damping  involved,  flutter  will 
occur  as  a net  result  of  energy  being  added  to  the  system.  However,  in 
reporting  his  results,  Whitehead  found  that  he  had  significant  mechani- 
cal damping  and  had  to  allow  for  an  average  value  as  a flutter  limit. 
This  limit  is  shown  in  Fig.  13.  If  the  value  of  the  imaginary  part  of 
the  moment  coefficient  is  below  this  lino,  then  flutter  is  predicted. 

I The  point  where  it  is  just  possible  for  the  flutter  to  exist  is  referred 

* to  as  a flutter  boundary.  Whitehead  observed  that  as  the  Mach  number  is 

^ increased,  the  points  at  which  flutter  is  just  possible  move  to  pro- 

1 gressively  lower  values  of  the  frequency  parameter.  This  corresponds  to 

I higher  fluid  velocities  or  to  lower  blade  stiffness.  The  effect  of  Mach 

number  is  therefore  highly  significant. 


( 


Figure  14  shows  the  frequency  parameter  below  which  torsional 
flutter  is  just  possible  and  has  been  plotted  against  blade  axis 
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Fig.  13  Torsional  Flutter  Boundaries  . 
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Fig.  14  Sub-Critical  Flutter 
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position.  Once  again,  a comparison  is  made  by  utilizing  the  Jones/Rao 
technique  and  the  Whitehead/Smith  technique.  The  results  are  in  close 
agreement  and  concur  with  Whitehead's  observation  that  increasing  the 
Mach  number  is  favorable  and  that  the  worst  position  for  the  torsional 
axis  is  around  fifty  to  seventy  percent  chord. 

Two  Degrees  of  Freedom 

A numerical  investigation  for  two  degrees  of  freedom  was  conducted 
using  the  program  for  unsteady  airload  prediction  and  interfacing  it 
with  the  program  developed  for  flutter  analysis.  The  geometrical  and 
structural  properties  of  the  cascade  reference  airfoil  used  in  this 
analysis  were  calculated  form  data  obtained  from  Refs.  18,  19,  and  20. 
These  properties  are  given  in  App.  B.  To  obtain  the  flutter  speed  and 
frequency,  the  program  was  executed  for  the  given  cascade  according  to 
the  procedure  as  outlined  in  the  block  diagram  of  Fig.  B-1  in  App.  B. 

Figure  15  represents  the  results  obtained  for  the  flutter  analysis 

for  various  values  of  interblade  spacing  (s  = l.Ot,  l.l?,  1.2)i,  and 

1.3£)  while  the  interblade  stagger  angle  was  varied  from  44  degrees  to 

60  degrees.  Over  the  range  of  values  considered,  it  can  be  seen  from 

Fig.  15  that  as  the  interblade  stagger  angle  is  increased,  the  flutter 

Mach  number  also  increases.  This  indicates  that  an  increase  in  the 

interblade  stagger  angle  would  be  beneficial  in  preventing  flutter  at 

the  lower  Mach  numbers.  This  trend  is  in  general  agreement  with  a re- 

21 

port  by  White  and  Rao  in  which  a cascade  of  helicopter  rotor  blades 
was  studied  at  various  stagger  angles  and  rotor  blades  spacings. 


Figure  16  shows  the  effect  of  varying  the  interblade  spacing  for 
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various  interblade  stagger  angles  (x  = 46^^,  52°,  56°,  and  60°).  The 
^nd  of  the  curves  shows  that  as  the  interblade  spacing  increases  for  a 
given  interblade  stagger  angle,  the  flutter  Mach  speed  decreases.  This 
indicates  that  increasing  the  interblade  spacing  will  have  a destabiliz- 
ing effect  on  the  cascade  of  blades.  This  same  trend  was  also  observed 
by  White  and  Rao  in  their  work  with  helicopter  rotor  blades. 

A summary  of  the  results  of  the  two  degree  of  freedom  flutter 
analysis  is  presented  in  Table  3. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

A unique  velocity  potential  lifting  surface  technique  has  been  used 
to  determine  the  unsteady  airloads  for  an  infinite  cascade  of  staggered 
blades  in  subsonic  flow.  The  unsteady  airloads  were  then  utilized  in  a 
FORTRAN  program  to  perform  a flutter  analysis  for  a single  degree  of 
freedom  system  in  torsion  and  for  a two  degree  of  freedom  system  in 
bending  and  torsion. 

Several  geometric  and  flow  parameters  of  a staggered  cascade  were 
varied  over  a specific  range  of  values  and  were  found  to  influence  the 
cascade  in  the  following  manner: 

1.  If  Mach  number  is  increased,  it  will  have  an  increased  damping 
effect  on  the  stability  of  the  cascade  in  pitching  motions. 

2.  If  interblade  spacing  is  increased,  the  damping  decreases  while 
still  below  the  first  critical  frequency  and  is  equal  to  zero 
at  the  first  critical  frequency. 

3.  As  interblade  stagger  angle  increases,  the  lowest  critical 
frequency  decreases. 

4.  As  interblade  phase  lag  is  increased  with  increasing  reduced 
frequency,  the  damping  will  decrease  rapidly  to  zero  for  an 
interbladc  phase  lag  of  180  degrees. 

When  the  cascade  was  analyzed  for  a single  degree  of  freedom  sy.tem 
in  torsion,  results  were  found  to  compare  favorably  with  the  results 
obtained  by  D.S.  Whitehead  for  flutter  boundaries.  A comparison  of  the 

B 

lifting  surface  technique  employed  by  Rao  and  Jones  and  a different 


I 


technique  utilized  by  Whitehead  in  determining  torsional  flutter,  were 
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shown  to  be  in  agreement  in  indicating  that  increasing  Mach  number  is  a 
desirable  condition. 

16  21 

By  utilizing  work  performed  by  Rao  and  White  and  Rao  , a FORTRAN 
' program  was  developed  that  is  capable  of  determining  the  flutter  speed 

and  frequency  of  a staggered  cascade  of  blades  in  unsteady  subsonic  flow. 
This  program  is  valid  for  a two  degree  of  freedom  system  in  bending  and 
torsion.  It  was  determined  that  increasing  interblade  stagger  angle 
would  be  beneficial  in  preventing  flutter  at  low  Mach  numbers  while  an 
increase  in  interblade  spacing  would  be  a destabilizing  influence  on  the 
cascade  of  blades. 

While  a new  flutter  program  was  developed  to  determine  the  cascade 
flutter  speed  and  frequency,  only  a limited  range  of  values  were  analyzed 
when  varying  cascade  parameters.  This  was  due  primarily  to  the  iterative 
scheme  that  was  involved  and  the  associated  expense  of  computer  compu- 
tational time.  Even  though  cost  effictiveness  was  considered  as  a criter- 
ion in  the  development  of  this  program,  the  program  still  remains  as 
being  a reasonably  efficient  technique  to  be  used  in  determining  the  un- 
steady  aerodynamic  derivatives,  the  flutter  speed,  and  the  flutter  fre- 
quency for  a given  cascade.  By  way  of  example,  the  run  time  associated 
^ with  determining  the  flutter  speed  and  frequency  for  the  example  given  in 

I App.  B,  utilized  less  than  twelve  seconds  of  CPU  execution  time  on  the 

AMDAHL  470V/6  computer  at  Texas  A^iM  University.  However,  further  improve- 
i ments  could  be  made  in  the  overall  efficiency  of  the  program  and  the  iter- 

ative procedure  it  utilizes.  Additionally,  investigation  should  also  be 
conducted  over  a larger  range  of  cascade  parameter  values  to  determine 
their  effect  on  the  cascade  flutter  speeds  and  frequencies. 
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APPENDIX  A 


Physical  Interpretation  of  Resonance  Conditions 


Consider  the  airfoil  as  (jiven  below  in  Fig.  A-l  , 


s cosX 


Reference  Airfoil 


s sin>. 


Fig.  A-1  Resonance  Model  . 

where  X is  stagger  angle,  s is  the  distance  between  the  leading  edge  of 
consecutive  airfoils,  and  U is  freestream  velocity.  The  influence 
created  by  the  unsteady  motion  of  the  airfoil  will  propagate  and  influ- 
ence the  reference  airfoil  after  time  t.  This  gives, 

[U  t]2  = [s(sinX)  - Ut]"  + [s(cosx)]2  (A1 ) 


U^t^  = _ 2Uts(sinX)  + U^t^  + s^{cos^x) 


[U^  - U^Jt^  + [2Lls(sinX)]  t - s^  = 0 


solving  for  t gives. 


(U2  - U2) 


t = - MsinA).. 


(U?  - U2) 
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given  that  6 = [l  - and  = - , Eq.  (M)  becomes, 


t = - -^-^(sinA)  1 /sin^A  + (-  - iT. 


From  Ref.  7,  it  is  shown  that, 

pt  + 0 = 0,  27r,  etc.  , (AG) 

where  P(=^)  is  the  frequency  of  blade  oscillation  in  rad/sec,  and  o is 
phase  lag.  Now  Eq.  (A5)  becomes, 

-o  = - — ajS(sinx)  + — - wS  /sin^x  + (-  - 1)  (A7) 

B^  V m2 

where  S(=|-)  is  space  to  blade  semi-chord  ratio.  After  consol i dati ng 
terms,  Eq.  (A7)  becomes. 


= ± /-  - cos^A  - sinA 

° V m2 


The  critical  or  resonance  frequency  is  then  given  as. 


w.  = ± /-  - cos2a  - sinX  ^ 

c V m2  ^ 
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APPiNDlX  B 

Example  Problem  for  Two  Degree  of  Freedom  Flutter 

A numerical  test  case  for  cascaded  airfoil  flutter  analysis  was 

coii;puted  using  a flutter  analysis  FORTRAN  program  interfaced  with  the 

unsteady  airload  prediction  program.  The  geometrical  and  structural 

properties  of  the  reference  airfoil  used  in  this  analysis  are  calculated 

from  data  obtained  from  Refs.  18,  19,  and  20,  and  are  given  below. 

Static  moment  about  ^ 

the  elastic  axis  S = 5.35x10”^  slug-ft. 

Moment  of  inertia  ^ 

about  elastic  axis  I = 6.316x10  ^ slug-ft^. 

Mass  per  unit  span 

of  blade  M = 0.0214  slugs/ft. 

Bending  natural 

frequency  = 1000  rad/sec. 

Torsional  natural 

frequency  = 2000  rad/sec. 

The  cascade  parameters  used  are. 

Blade  semi -chord 

length  i = 0.0833  ft. 

Interblade  spacing  s = 1 .2j!, 

Interblade  phase 

lag  6 = -0.0835 

Interblade  stagger 

angle  X = 52  deg. 

To  obtain  the  flutter  speed  and  frequency,  the  program  was  executed 
for  the  given  cascade  according  to  the  procedure  as  given  'n  the  block 
diagram  of  Fig.  B-1.  After  the  input  parameters  arc  read  into  the  pro- 
gram and  initial  values  are  assunied  for  the  Mach  number  and  the  reduced 
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Fig.  B-1  Program  Flowchart  . 
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frequency,  the  program  calculates  the  unsteady  aerodynamic  derivatives 
tiiat  are  to  be  used  as  the  forcing  functions  in  Lagrange's  equations  of 
motion.  The  aerodynamic  derivatives  that  are  initially  obtained  are 
referred  to  the  blade  mid-chord  axis.  Since  it  is  necessary  for  the 
aerodynamic  derivatives  *'0  be  referred  to  tlie  elastic  axis,  a coordinate 
transformation  is  performed  as  explained  in  App.  C.  After  the  coordinate 
transformation,  the  flutter  determinant  is  solved.  This  yields  PI  and  P2 
cs  solutions  to  the  quadratic  equation  of  tiie  real  part  of  the  flutter 
determinant  and  P3  as  a solution  to  the  linear  equation  of  the  imaginary 
part  of  the  flutter  determinant.  Since  an  iterative  scheme  is  involved 
in  obtaining  the  flutter  speed  and  frequency,  initial  values  of  Mach 
number  and  reduced  frequency  are  assumed.  The  initial  values  of  Mach 
number  are  0.500  and  0.5100.  Several  values  of  w are  then  calculated  for 
each  Mach  number  and  tfie  respective  values  of  PI,  PZ , and  P3  are  found. 
These  values  are  then  plotted  as  shown  in  Fig.  B-2.  Figure  B-2  indicates 

that  P2  is  equal  to  P3  when  (u  = 0.3600  for  M = 0.500  and  Fig.  B-3  shows 

that  P2  is  equal  to  P3  when  u = 0.3486  for  M = 0.510.  These  values 

represent  possible  flutter  frequencies  if  the  assu'iied  .^lach  number  is 

equal  to  the  flutter  Mach  number.  When  the  assumed  Mach  number  is  equal 
to  0.500,  Fig.  B-2  shows  02  = P3  = 2293  at  w = 0.3600.  If  these  values 
are  used  in  the  following  equation, 

K = (Bl) 

ua 

where  a = 1037.26  ft/sec.  , then  the  flutter  Mach  number  should  be 
equal  to  0.5116.  However,  since  this  is  not  equal  to  the  assumed  Mach 
number  of  0.500,  another  iteration  must  be  performed.  The  next  ite-ation 
for  an  assumed  Mach  number  of  0.510  is  shown  in  Fig.  B-3.  Here  it  can 
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be  seen  that  P2  = P3  = 2290  at  w = 0.3486.  Again,  using  Eq.  (Bl),  the 
flutter  Mach  number  is  calculated  to  be  0.527G.  Once  again,  the  calcu- 
lated Mach  number  is  not  equal  to  the  assumed  Mach  riumber,  so  a flutter 
speed  still  has  not  been  determined.  Tor  the  next  Mach  iteration.  Fig. 
B-4  will  provide  the  assumed  Mach  number.  If  the  p;revious  assunied  Mach 
numbers  are  plotted  versus  the  calculated  flutter  Mach  numliers  , it  can 
be  seen  from  Fig.  B-4  that  an  extrapolated  value  of  M = 0.480  should  be 
used  for  the  next  assumed  Mach  iteration.  If  this  is  done,  then  P2'-r3  = 
2297  at  u = 0.3837.  Using  Eq.  (Bl),  a calculated  flutter  Mach  number 
value  of  M = 0.4809  is  obtained.  Since  this  value  is  essentially  the 
same  as  the  assumed  Mach  number,  it  is  therefore  considered  a valid 
flutter  apeed  at  a flutter  frequency  of  = 0.3837. 


APPENDIX  C 

Change  of  Reference  Axis 

12 

Jones  has  given  a convenient  axis  transformation.  In  Fig.  C-1 , 

I 

let  the  reference  axis  be  moved  from  mid-chord  to  a point  0 , at  a dis- 
tance hJc,  forward  of  mid-chord. 


o'  0 


z z 


Fig.  C-1  Blade  Reference  Axis. 

When  referred  to  the  original  coordinates  z and  a,  the  displacement 
of  a point,  p,  downwards  is, 

Z = P.Z  + tCa  + «,(z  + fa)  (Cl) 


(C2) 

(C3a) 

(C3b) 
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If  L and  M are  the  lift  and  moment  referred  to  0 (rni d-cfiord)  and 

I 

L amd  M are  the  corresponding  force  and  moment  referred  to  U , i t can  be 
shown  that  the  relation  exists  as  shown  in  Fig.  C-2. 


L 


hi  — 


M 


Fig.  C-2  Force  and  Moment  Orientation  . 

From  Fi g.  C-2 , 

r = L (C4a) 

M = M - h?.L  (C4b) 

where  M and  M are  positive  in  the  clockwise  direction.  It  can  be  shown 
that. 


_L_  = _L  _ 

pU^Ji,  pU^Jl 


(C^zR 


'tzi 


triR 


'f.al' 


Si  mi  larly. 


^^izR  ^ ^ • 


" ^^mzR  ^ ^ ^ • 


(C5a) 

(C5b) 

(C5c) 


ice) 


From  Fqs.  (C5c)  and  (C6),  the  expressions  for  the  new  aerodynamic 
derivatives  are  given  in  terms  of  the  original  ones  referred  to  the 
mid-chord.  These  equations  are  given  as. 


^i:zl  " ^iZl 


''taR  ' ^VaR  ' ^^fczR 


'^iul  ^Jal  ' ^^f  Zl 


^iiizR  ' ^inzR  ^ *^^tzR 


^mzl  ^mzl  ^^fzl 


'^nuR  ' ^mcR  ' ^^^mzR 


^z.r)  - 


tzR 


- h(C_i  - C^,.,)  - h^C, 
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(C7a) 

(C7bj 

(C7c) 

(C7d) 

(C/e) 

(C7f) 

(C7g) 

(C7n) 


^nul  ^nul  "'"nizi  ^ ''tal'  ' " ^'£zl 
The  equations  in  (C7)  are  a linear  conin’nation  of  the  original 
derivatives  referred  to  the  mid-chord.  These  equations  apply  '/iien  the 
reference  axis  is  moved  aft  of  the  mid-chord  a distance  ht  towards  the 
trailing  edge. 
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